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Abstract 



We consider the ground state in a model with scalar nucleons and a meson 
using the formalism of light-front dynamics. Light-front potentials for two- 
nucleon bound states are calculated using two approaches. First, light-front 
time-ordered perturbation theory is used to calculate one- and two-meson- 
exchange potentials. These potentials give results that agree well with the 
ladder and ladder plus crossed box Bethe-Salpeter spectra. Secondly, approx- 
imations that incorporate non-perturbative physics are used to calculate al- 
ternative one-meson-exchange potentials. These non-perturbative potentials 
give better agreement with the spectra of the full non-perturbative ground- 
state calculation than the perturbative potentials. For lightly-bound states, 
all of the approaches appear to agree with each other. 

PACS number(s): 21.45. -Fv, 03.65. Ge, 03.65.Pm, ll.lO.Ef 
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I. INTRODUCTION 



Recent experiments at Thomas Jefferson National Accelerator Facility have measured 
the A(Q^) structure function of the deuteron for momentum transfers up to 6 (GeV/c)^ [|l|, 
and measurements for B{Q^) are planned. At such large momentum transfers, a relativistic 
description of the deuteron is required. One approach that gives such a description is 
light-front dynamics, which we will examine here. To separate the effects of the using light- 
front dynamics from the effects of the model, we choose to use the massive Wick-Cutkosky 
model. This is a "toy model" investigation, instead of the full nuclear theory calculation. 
Using this model, the light-front Hamiltonian approach is used to solve for the bound-state 
wavefunction. The results of our calculation can then be compared to other calculations 
done with the same model but different approaches. The simplest observable that can be 
compared is the relationship between the bound-state mass and the coupling constant. 

The utility of the light-front dynamics was first discussed by Dirac |^. We start by 
expressing the four-vector in terms of the light-front variables 

where = ± x^. This is simply a change of variables, but an especially convenient one. 
Using this coordinate system and defining the commutation relations at equal light-front 
time {x~^ = tLp), we obtain a light-front Hamiltonian The Hamiltonian is used in the 

light-front Schrodinger equation to solve for the ground state. 

There are many desirable features of the light-front dynamics and the use of light-front 
coordinates. First of all, high-energy experiments are naturally described using light-front 
coordinates. The wave front of a beam of high-energy particles traveling in the (negative) 
three-direction is defined by a surface where x^ is (approximately) constant. Such a beam 
can probe the wavefunction of a target described in terms of light-front variables 0,^: the 
Bjorken x variable used to describe high-energy experiments is simply the ratio of the plus 
momentum of the struck constituent particle to the total plus momentum (p"^) of the bound 
state. Secondly, the vacuum for a theory with massive particles can be very simple on 
the light front. This is because all massive particles and anti-particles have positive plus 
momentum, and the total plus momentum is a conserved quantity. Thus, the naive vacuum 
(with p+ = 0) is empty, and diagrams that couple to this vacuum are zero. This greatly 
reduces the number of non-trivial light-front time-ordered diagrams. Thirdly, the generators 
of boosts in the one, two, and plus directions are kinematical, meaning they are independent 
of the interaction. Thus, even when the Hamiltonian is truncated, the wavefunctions will 
transform correctly under boosts. Thus, light-front dynamics is useful for describing form 
factors at high momentum transfers. A drawback of the light-front formalism is that the 
Hamiltonian is not manifestly rotationally invariant, since the generators of rotations about 
the one and two directions are dynamical. A study of the effects of the loss of rotational 
invariance of the excited states in the model being used here was made in Ref. [^, which 
shows that there is less breaking of the degeneracy in the spectrum when higher order 
potentials are used as opposed to lower order potentials. 

There are other approaches that can be used to obtain relativistic wavefunctions and 
bound-state energies, including the Feynman-Schwinger representation (FSR) of the two- 
particle Green's function and the Bethe-Salpeter equation (BSE) [^-[1^. The FSR is 
useful since it can be constructed so it is equivalent to the Bethe-Salpeter equation using a 
kernel where a// two-particle-to- two-particle ladder and crossed ladder diagrams are included. 
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However, it requires a path integral to be done numerically, so it is computationally intensive 
to obtain an accurate answer. The FSR result is to be considered the full solution that the 
Bethe-Salpeter and Hamiltonian equations approximate. The BSE can be solved much 
quicker than the FSR, however a truncation of the BSE kernel is required, causing the BSE 
results to differ from the FSR results. It is well known that any finite truncation of the 
kernel yields bound-state wavefunctions with problems, such as the incorrect one-body limit 

. For our scalar model the potential is always attractive, so truncation of the kernel gives 
binding energies that are too small |1T^-|T6|. Another approach is explicitly covariant light- 
front dynamics |]T7|,|T^ , where manifest covariance is kept at the price of using a null-plane 
whose orientation is not fixed. 

Here we study light-front dynamics because of its close connection to experimental ob- 
servables. Using field theory, light-front potentials can be derived that give results physically 
equivalent to those of the Bethe-Salpeter equation. We define "physically equivalent" in sec- 
tion |111 (J|. Depending on the diagrams used to construct the potential, one can argue that 



certain potentials are physically equivalent to the Bethe-Salpeter equations with certain ker- 
nels. The best that these potentials can do is reproduce the results of the corresponding 
Bethe-Salpeter equations. On the other hand, light-front potentials can be constructed that 
attempt to incorporate non-perturbative physics. It is possible that these potentials can 
give results that agree with the full theory results better than the Bethe-Salpeter equation 
using low-order kernels. We will consider both types of potentials in this paper, and see how 
well they perform. 

A brief discussion about the approximation used is in order. It is well known that the 
vacuum of the full Wick-Cutkosky model is unstable [|I^ due to the cubic coupling which 
provides the interaction. However, when the bound-state calculation is restricted to the 
two-particle sector, the quenched approximation is used, and the self-energy and vertex- 
renormalization diagrams are neglected, the theory has a well defined ground state. In 
this paper, we compare the results of our light-front Hamiltonian calculation to the Bethe- 
Salpeter and FSR calculations, both of which use the same approximations. The use of 
these simplifying approximations allows us to highlight the differences between the various 
approaches. The inclusion of the self-energy diagrams and counterterms for the light-front 



Hamiltonian 2G and for the FSR 21,22 will not be discussed here 



A. Outline of the paper 

The objective of this paper is to obtain the bound-state energy for the ground state 
in our theory. In our model, neutral scalar nucleons interact via a Yukawa interaction, 
which is mediated by a neutral scalar meson. The light-front Hamiltonian derived from the 
Lagrangian is used in the light-front Schrodinger equation for a two-nucleon bound state. 
The rules for light-front time-ordered perturbation theory (LFTOPT) are then derived for 
this Schrodinger equation, along with the Feynman rules for the effective potential. All of 
this is discussed in more detail in section y. 

The full potential for the light-front Schrodinger equation is given by an infinite sum 
of diagrams. Using the LFTOPT rules in section |ITT|, we derive the one-boson-exchange 
(OBE) and two-boson-exchange (TBE) potentials, where the diagrams that give rise to 
mass and vertex renormalization are not included. The two-boson-exchange stretched-box 
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(TBE:SB) diagrams, a subset of the TBE diagrams, are used to construct the TBE:SB 
potential, and the utihty of this potential is commented on. After the discussion of the per- 
turbative potentials, we define several potentials in section |rV| which attempt to incorporate 
non-perturbative physics in a OBE type of potential. Three potentials are obtained by ap- 
proximating the OBE potential directly (giving the symmetrized mass, instantaneous, and 
retarded potentials), and one potential is obtained by approximating the Bethe-Salpeter 
equation and reducing it to a three-dimensional equation (giving the mo dified- Green 's- 
function potential). 

In section |V], the potentials are used to numerically obtain the spectra, the coupling con- 
stant versus bound-state mass curves. This is done by solving the light-front Schrodinger 
equation with each of the truncated potentials (OBE, TBE, and TBE: SB), and the ap- 
proximate potentials (symmetrized mass, instantaneous, retarded, and mo dified- Green 's- 
function). We compare these results to those in the literature obtained with other ap- 
proaches. 

We summarize our findings in section Only a few low-order terms from the pertur- 
bative potentials are needed to approximate the results of the physically equivalent Bethe- 
Salpeter equations [ff|, p3| , p4| . Since the light-front potentials are calculated without including 



the mass and vertex renormalization diagrams, the Bethe-Salpeter kernels used also do not 
include those diagrams. Since the interaction in the Wick-Cutkosky model is strictly at- 
tractive, the spectra calculated using the perturbative potentials will underestimate the the 
binding energy compared with the spectra for the physically equivalent Bethe-Salpeter equa- 
tions. As progressively higher-order terms in the potential are kept, the spectra calculated 
will agree better with the BSE spectra. However, when a truncated kernel is used in the 
Bethe-Salpeter equation, the solutions obtained are known to be a poor approximation of 
the full solution which for our theory (where no renormalization graphs are kept) is 



given by the Feynman-Schwinger representation of the Green's function ||8|j25[]. Hence, our 
truncated potentials cannot give the full results. The non-perturbative approximations of 
the potential give spectra that more closely match the spectrum of the full solution. This 
suggests that these approximations are reproducing the physics more accurately than the 
perturbative potentials. 

The conventions and notations employed in this paper are summarized in Appendix 
^ The conversion of the light-front Schrodinger equation into matrix form suitable for 
numerical evaluation is discussed in Appendix |B[ Azimuthal-angle integrations of the OBE 
potentials, which help simplify the evaluation of the bound-state wavef unctions, are given 
in Appendix 0. The loop integrations and azimuthal-angle integrations needed for the 
TBE potentials are discussed in Appendix 0. A check of the validity of the uncrossed 
approximation used in section [IV D| is done in Appendix |^. 

This study is related to the work of Sales et al. ||23|, who computed bound states using 
the light-front Hamiltonian with the OBE and TBE:SB approximations, and compared to 
the ladder Bethe-Salpeter equation results. Here, we consider also the effect of including the 
crossed graph part of the TBE potential as well as several non-perturbative approximations. 
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II. OUR MODEL 



We consider an isospin doublet of two uncharged scalars (f) = (0i,02) with mass M 
(which we will refer to as nucleons), that couple to a third, uncharged scalar x with mass n 
(which we will refer to as a meson) by a (p'^x interaction. This is the massive extension of 
the Wick-Cutkosky model [p6|, which has been used on the light front to study scattering 



states as well as bound states |23]. The Lagrangian is 



(1) 



where 5^ is a dimensionless coupling constant and 0^ = 0^ + 02- 



A. Light-front Hamiltonian 



To obtain the light-front Hamiltonian from the Lagrangian in Eq. (|l]), we follow the 
approach used by Miller P] and many others (see the review f^) to write the light-front 
Hamiltonian (P~) as the sum of a free, non-interacting part and a term containing the 
interactions. We use the conventions given in Appendix 0. The operators we use can be 
expressed in terms of Fock space operators since for this theory in light-front dynamics, the 
physical vacuum is the Fock space vacuum, and thus the Hilbert space is simply the Fock 
space. The Hamiltonian is obtained by using the energy-momentum tensor in 



The usual relations determine T^^, with 



dC 



d{d^(t)r) 



(2) 



(3) 



in which the degrees of freedom (the fields and x) ^^re labeled by ipr- 

It is worthwhile to consider the limit in which the interactions between the fields are 
removed. This will allow us to define the free Hamiltonian Pq and to display the necessary 
commutation relations. The energy-momentum tensor of the non-interacting fields is defined 
as Tq" . Use of Eq. (|^) leads to the result 



rp^U 



with 



d''(j)d''(f) - 



-t n 
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-M 



2 J,2 



+ d^xd^X 
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+ M 



2 J,2 



2,, 2 



V_lX ■ + fJ' X 



(4) 



(5) 



The scalar nucleon fields can be expressed in terms of creation and destruction operators: 



(27r)3/2V2Ap 



ik-x 



(6) 
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M2 + fc2 



where z = 1,2 is a particle index, k ■ x = ^{k~x~^ + k^x~) — kj^ ■ x± with k~ — k+ i 
and k = {k~^, k±). Note that k~ is such that the particles are on the mass shell, which is a 
consequence of using a Hamiltonian theory. The 6 function restricts k~^ to positive values. 
Likewise, the scalar meson field is given by 



(Pk^dk^ e{k+) 
(27r)3/2v/2A^ 



a^{k)e~'''"' + al{k)e 



Ak-x 



(7) 



where k~ 



^+ so that the mesons are also on the mass shell. The non- vanishing 
commutation relations are 



(fc),al(fc')l =6{k^-k'^_)5{k+ -k'^ 



(8) 



where a = 1, 2, x is a particle index. The commutation relations are defined at equal light- 
front time, x"*" = 0. It is useful to define 



- k') = S{k^ - k'^)5{k+ - k'+), 

which will be used throughout this paper. 

We write a ket in the two-distinguishable-particle sector of the Fock space as 



(9) 



\kuk2) = a\{ki)al{k2)\0). (10) 
This implies that the identity operator in this Fock space sector can be written as 

l2 = J d^ki^^dkf J d'^k2,±dk^\ki,k2){ki,k2\. (11) 

The derivatives appearing in the quantity Tg^^ are evaluated and then one sets x~^ to 
to obtain the result 



M2 + fci / ^ 



a\{k)ai{k) + a2{k)a2{k)j + 



;/2 4- 

i^^a{{k)a,ik) 



(12) 



with Jj. = J d'^k±dk~^ 6{k^). Eq. (|T2D has the interpretation of an operator that counts the 



for the nucleons and 



k+ 



for the mesons) of all of 



light-front energy k~ (which is — ^ 
the particles. 

We now consider the interacting part of the Lagrangian, Cj. An analysis similar to that 
for the non-interacting parts yields the interacting part of the light-front Hamiltonian Pf; 



^ 2 JkJk' {2n)y^^2k+k'+{k+ + k'+) 



i=l,2 



X { [2al{k + k')a^{k')ai{k) + a^^{k + k')ai{k')ai{k) 
+ Hermitian conjugatej. 



(13) 



The interaction Hamiltonian is self-adjoint since the Hilbert space is the Fock space. The 
total light-front Hamiltonian is given by P~ = Pq + gPf. 
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B. Hamiltonian bound-state equations 



We will be studying the bound states of two distinguishable nucleons. The technology of 
time-ordered (old-fashioned) perturbation theory is used to construct the light-front time- 
ordered perturbation theory (LFTOPT) for our Hamiltonian. We start with the light-front 
Schrodinger equation in the full Fock space, 

{p,-+gPf)\^3^) = \^P^^)P^s. (14) 

where Pg" + d^F is the Hamiltonian in the full Fock-space basis, \ipp^) is the ground-state 
wavefunction in the full Fock space, and P^g is the light-front energy of that state. Recall 
that Pq , the non-interacting part of the Hamiltonian, is diagonal in the momentum basis, 
while Pf , which contains the interaction, has only off-diagonal elements. 

A serious drawback of this equation is that the wavefunction ipp^) has support from 
infinitely many sectors of the Fock space, since Pf changes the total number of particles. 
However, the components of the wavefunction with many particles will be small compared 
to the two-particle component if the coupling constant is not too large. We will construct 
the two-particle light-front Schrodinger equation which the two-particle component of the 
wavefunction satisfies. From this construction, we will obtain the rules for the LFTOPT. 

We start by introducing the projection operators V and Q. The operator V projects out 
the sector of Fock space with two distinguishable nucleons and no mesons, while Q = I — V 
projects out all the other sectors. We define 

V\i/jf) = IV^^S) (15) 

Qli^r) ^ l^g'), (16) 

so that {ipW) — 14''^^) + I'^Q^)- Since the free Hamiltonian does not change the number of 
particles, [V, Pq] = [Q, Pq] = 0. The interaction Hamiltonian changes the particle number, 
so it cannot connect the two-particle sector to itself, thus VPfV = 0. 

Using these projection operators, Eq. (IMD can be broken up into two parts. 



Pol^""') + gVPfQl^^') = |^^')Pgs (17) 
Pq + gQPFQ) li'f) + gQPfVli:''^) = |^g^)PGs- (18) 



Eliminating the \ipQ^) and using the expression of the identity given in Eq. ([TT]) we obtain 
the two-particle effective light-front Schrodinger equation 

J (fpi,^dpt J d^P2,Ldpt{ki,k2\ [Pq + V{g,PGs)] \Pi,P2){pi,P2\^'^^) 

= (fcl,fe2|^^')PGS, (19) 

where Pq and the potential V act in the two-nucleon basis. The two-particle potential is 
given by 

v(,. P-) = ,^vprp—p^^j^prv. (20) 
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Note that Eq. ([T9|) is similar to Eq. (0), except for two main differences. Here we have 
a two-nucleon wavefunction, which makes it simpler. However, the potential is light-front 
energy dependent, which makes it more complicated. 

The denominator in the definition of the potential is non-diagonal in the full Fock space, 
so the matrix inversion that it represents is highly non-trivial. This problem is avoided by 
expanding the inversion in powers of the coupling constant g to get 



v{g,p-) = vpr 



g'Q g (p- gQ 



PfV. (21) 



This can be simplified further by noting that in the two-nucleon sector of our theory, every 
meson emitted must be absorbed, so there must be an even number of interactions. Thus, 
the full potential can be written as the sum of n meson exchange potentials, 

oo 

ViP-,g) = Y.9'''V^2n)iP-), (22) 

n=l 

where V(2n) is the potential due to the exchange of n mesons, given by 

V(2n)(P-)=P(^Pf ^__^_ j PfV. (23) 

It is easy to see how to write a sum of diagrams for the potential when one says what 
Eq. (|23| ) represents in words. We start off with two particles, then the interaction occurs. 
There are two possibilities of what can happen; nucleon 1 or 2 can emit a meson. Each 
possibility has a separate diagram. After the interaction, there is propagation with the 
light-front Green's function, 

GLFiP-) = p-\p- ^ (24) 

until another interaction occurs, and so on. We simply sum up all of the possible orderings 
of the interaction to get the full potential. The v}^ order potential is simply the sum of all 
possible diagrams with n-meson exchanges. 

Each intermediate state in Eq. (|23| ) has more than two particles, so the diagrams are 
two-particle irreducible with respect to the two-particle Green's function G'2LF = VGi^pV. 
We can represent G'2LF by its diagonal matrix elements, 

G2lf(/ci, fcs; P-) = p_ ,^ (25) 

In the diagrams we draw, the nucleons will be represented by solid lines and the mesons 
by the dashed lines. Although the states we will be considering consist of two distinguishable 
nucleons, we will not label the nucleon lines. We will be using the quenched approximation 
(so there are no nucleon loops) and neglect the mass and vertex renormalization diagrams 
(so the physical masses and coupling constant are used, and each meson emitted from one 
nucleon must be absorbed by the other nucleon). It is not expected that these restrictions 
will lead to qualitatively different results than the true full solution when the states are 
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not too deeply bound. The quenched approximation is reasonable when the masses of the 
nucleon fields are large compared to the binding energy. Use of the physical masses and 
coupling constant are reasonable as well when the momenta are not too large. 

We stress again that we will compare various truncations of the light-front Hamiltonian 
to other calculations which do not include renormalization diagrams. This is because we 
want to determine the effect of truncation on the light-front Hamiltonian. The differences 
between our calculation and those which include the self-energy graphs which may be 
large for deeply-bound states, are not considered here. 

The rules for drawing the n-meson-exchange graphs that correspond to this approxima- 
tion are: 

1. Draw all topologically distinct time-ordered diagrams with n mesons. Use solid lines 
for the nucleons and dashed lines for the mesons. 

2. Delete all graphs which couple particles to the vacuum. In the massive theory we 
consider here, these diagrams always vanish since the vacuum has zero plus momentum, 
and massive particles always have positive plus momentum. 

3. Our quenched approximation and use of the physical masses and coupling constant 
requires us to delete all graphs that have nucleon loops or have mesons that are emitted 
and absorbed from the same nucleon. 

4. Delete all other graphs which are not allowed in the particular approximation that 
is being considered. For example, consider the potential from the the Hamiltonian 
theory that can be obtained from the ladder Bethe-Salpeter equation. That potential 
will not have any graphs where the meson lines cross. 

Once the diagrams are drawn, we use the following rules to convert the sum of diagrams 
into the potential (fci, fc2|V(2n)(-P^)|Pi,P2): 

1. Overall factor of '"'"^(fei+fe2-pi-p2) ^ This delta function says that the total light-front 

2(2^)3 ^fc+fe+p+p+ 

three-momentum is conserved. We define the light-front three-momentum P = ki+k2. 

2. To each internal line, assign a light-front three-momentum Qi where i = 1,2, ... ,N 
and is the number of internal lines. The light-front energy for particle i with mass 

rrii is Qi = It is useful to define Zi = qf / . 

3. A factor of ^ for each internal line. 

z 

z 

4. An extra factor of p+^p- for each internal meson line. 

5. A factor of / — ry between consecutive vertices, where the sum is over only the 
particles that exist in the intermediate time between those vertices. 

6. Use light- front three-momentum conservation to eliminate all the independent mo- 
menta. 
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7. Integrate with / 2P+p-{2n)'i remaining free internal momenta. 

8. Symmetry factor of | when two nucleons are created or destroyed at the same time. 
With these rules, one can calculate the effective potential for any order. 

C. Further development of the Ught-front Schrodinger equation 

Once the potential is calculated, we can plug it into Eq. (p!9|), which we as 

/ d^Pi,±dpt J d^P2,±dpt{k^,k2\ [Po +V{g{p-),P-)] |pi,P2)(Pi,P2|^''') 

= (fci,fc2|V^^')P", (26) 

where P~ is an arbitrary light-front energy and g{P~) is the coupling constant which yields 
the bound-state wavefunction with P~ as the bound-state energy. We call this g{P~) the 
spectrum of the light-front Schrodinger equation for the corresponding wavefunction. 

The total momentum P = fci + ^2 is conserved by the potential given in Eq. (p2|), 
so the wavefunction in Eq. can be parameterized by the total momentum. To make 
the calculations easier later, we choose to be in the center-of-momentum frame, where the 
components of the total momentum can be written as P± = and P~^ = P~ = E. The 
ground-state energy, E, is the same as the mass of the bound state. In terms of the binding 
energy B, E = 2M — B. In the center-of-momentum frame, the ground-state wavefunction 
is parameterized by i?, so we can define 

(fci, fcsl^^/) = + fc2 - P)i^''\k,) (27) 

(fcl, fc2|V(2„)(P-)|pi,P2) = 5^''+^^! + k2-pi- P2)V^2n)iE; k,-p,). (28) 



With these, Eq. (p6|) effectively becomes a one-particle equation, where particle 2's momen- 
tum is determined by ^2 = P — ki. The minus component (the light-front energy) of particle 
2 is defined by the requirement that the particle 2 is on mass shell, so = {W^ + kl^^/k^. 
We also define x = kf / P"*" = xbj , where xbj is the Bjorken x variable, so that A;^ / P^ = 1—x. 
Likewise, we write the Bjorken variables that correspond to the momenta pi and qi as 
y = Pi / P^ and z = ql / P^ . 

Using Eqs. (P?]), (^), and the fact that the plus momentum of both nucleons is positive, 
we can write the light-front Schrodinger equation Eq. (pBD as 

/ d'p,,^£dptVigiE),E; k,; p,)^""^ (p,) = ^«S(fci)(P - k^ - k^). (29) 

It useful to convert from light-front coordinates ki = [k^ ,k±) to equal-time coordinates 
^ET = {k±, k^), using an implicit definition of k^ |27 



kf = ^ [fc°(fcET) + k^] (30) 



A;°(fcET) = + (31) 
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Often the explicit dependence of k'^ on A^et will not be shown. It is worth emphasizing 
that this is just a convenient change of variables; ip'^^lk-^T) is not the usual equal-time 
ground-state wavefunction. With this transformation, we can express kf and kf as 



kt 



1 + 



A;0 



E 
2"' 



(32a) 



''^ ~ (2^) 



Using these, Eq. (^91) becomes 



V{g{E), E, kET, Pet)V^'''(Pet) = ^''''(fcET) E'' - (2A;' 



GS/ 



(32b) 



Now consider the exchange of the particle labels 1 and 2. This causes 



k 



2,± 



kt 



= -ki,± 
E-kf 



(33) 



(34) 
(35) 



-k 



ET- 



which means that k^ as defined in Eq. (^) transforms as k^ —k^, so fcET - 
Consequently, exchange of particle labels 1 and 2 is the same as parity for ^et- 

Since the two nucleons are identical except for the particle label, the effective potential 
commutes with parity to all orders in g"^. Furthermore, the light-front Hamiltonian is explic- 
itly invariant under rotations about the three-axis. These considerations allow us to classify 
the wavefunctions as eigenfunctions of parity and the three-component of the angular mo- 
mentum operator. The ground state will have even parity and be invariant under rotation 
about the three-axis, so we can write 



^""""{kET) = i^^'^ikET, 4, 0fc) = i^^'^ikET, = ij^'^ikET, VT - 4), (36) 

where 6k and 0^ are the polar and azimuthal angles, respectively, for the vector A^et- 

The cylindrical symmetry and parity of the wavefunction can be used to rewrite Eq. (1331 

as 



dpET 
JO 



2ptptplT sin 



d9, 



V+ikET,ek;PET,Op)ij''''ipET,0. 



V''''(fcET,4) 



0\2 



E' - 4(A;") 



where 



{kET, Ok, Pet, Op) 

V{kET,Ok;PET,Op) 



V{kET, Ok;PET, Op) + V{kET, Ok, Pet, tt - Op] 

-V{kET;PET)- 



27r 



(37) 

(38) 
(39) 



We call V {kET, Ok', Pet, Op) the azimuthal-angle-averaged potential. The light-front vectors 
can also be expressed in terms of the azimuthal angle (p, where k = {k~^,kj_,(f)k), which 
allows the azimuthal-angle-averaged potential to be written in light-front coordinates 
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V{kt,h,r,Pt,Pi,±) = fr^^^V{k,-p,). (40) 

J Jo Zn 

All of the simplifications of Eq. (pTj) based on physical considerations have been ad- 
dressed. However, further rearrangements need to be done before Eq. (|37D is fit to be 
solved on the computer. Since these involve only numerical techniques, they are relegated 
to Appendix 0. 



III. PERTURBATIVE POTENTIALS 

A truncation must be made of the expansion of the potential given in Eq. (^21) , since it 
is not feasible to calculate the infinite sum of graphs for the potential in this Hamiltonian 
theory. In this paper we consider three truncations of the potential derived from the field 
theory. First the OBE potential and the TBE potentials are calculated. Then, we note that 
a subset of the TBE diagrams, the stretched-box diagrams, correspond to the truncated 
potential derived from the ladder Bethe-Salpeter equation. Thus, three truncated potentials 
are obtained that have a physical interpretation. 

The matrix elements of these potentials are written in the two-particle momentum basis, 
denoting the momentum of the incoming particles by pi and p2, and the outgoing particles ki 
and For simplicity, we choose to work in the center-of-momentum frame. By inspecting 
the rules for converting a light-front time-ordered diagram into a potential given in section 
[II B| , and looking at Eq. (|28|), we find that each piece of the effective-one-particle potential 
V{E; ki,pi) is proportional to 



2i27rr ^kfktptpt 

(An extra factor of E is included to simplify later equations.) This term will by suppressed 
in all of the potentials written in this paper. 



A. OBE Potential 



We start by drawing all the allowed and non-vanishing time-ordered diagrams with one 
meson exchange. These diagrams are shown in Fig. |I|. The light-front time-ordered per- 



turbation theory rules given in section [II B| are used to calculate the potential due to OBE 
potential. 



Vobe(^; ki;pi 



My 
'eJ 



9{x - y)/\x - y\ 



E -p^ -k^ - id (fci - pi) 
^ Ojy - x)/\y - x\ 
E -ki -p2 -uj-{pi - ki] 



(42) 



We have introduced the notation that meson with light-front three-momentum q has a 
light-front energy given by 
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VoBE 

ki - Pi / + \pi - ki 

—d ^ 

ki Pi ki Pi 

FIG. 1. The diagrams for the OBE potential. 



(43) 



The azimuthal-angle average of Vqbe is discussed in Appendix p. 

The potential given in Eq. ( ^21) can also be used for scattering states. In that case, 
E = ki + k2 = Pi + P2 , which allows the potential to be written as 



{ki -piY - /i^' 



(44) 



The scattering potential is the same as the usual equal-time OBE potential. This must be 
the case, since the scattering potential is also given by covariant Feynman diagrams, which 
have the same form independent of the form of dynamics. 

Returning to the bound-state regime, we note that the OBE potential can easily be 
written in terms of the equal-time coordinates. A reorganization of Eq. (^) yields 



VbBE(-E'; ^et; Pet) 



e) 



d{x - y) 



(kt - pt)iE -Pi- k^) - /i^ - (p^ - k^y 

e{y - x) 



(pt - kt) (E - ki -p-)-^2_ (p^ _ fc^)2 

Using the relations in Eq. (|32|) , we find 

'E^-AM^ ^ 



(45) 



{kt-pt){E-pi-k2 



' ^3 p3' 



Pet ^et 



(46) 



Under the exchange k ^ p, the only thing that changes in Eq. (^) is that the first term 
picks up a minus sign. This observation allows Eq. (|45|) to be rewritten as 



VoBE{E]kET;PET) 



E 



fc£ 



where 



E^ - AM^ _ 2 _ - 2 
9et ~ ^et ~ Pet- 



(47) 



(48) 
(49) 
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Note that g^^ is not the hght-front energy of the meson, since in a Hamiltonian theory only 
the hght-front three-momenta are conserved; the four-momenta is not conserved. Equation 
(^) will be useful in the context of approximations based on the physical arguments that 



we will discuss in section IV. 



B. TBE Potential 



As in the previous section, we start by drawing all the allowed, non-vanishing time- 
ordered diagrams with two meson exchanges shown in Fig. ^ The diagrams are classified 
according to the behavior of the intermediate particles. The total TBE potential is given 
by the sum of all the diagrams, so 



TBE 



TBE:SB 



TBE:SX 



TBE:TX 



TBE:WX 



TBEiZX- 



(50) 



In the diagrams for the TBE potential in Fig. |^, the intermediate loop momenta can be 
parameterized by qi or q2. The dependent variable is defined by the relation P = qi + q2. 
The Bjorken x variable that corresponds to qi (^2) is labeled with z {1 — z). We use the 
Feynman rules to calculate all of these potentials, starting with 



Vtbe:Sb(-E; ki;pi] 



_cPqj__ 
2(27r)= 



dz 



X ■ 



9{z-y)9{x-z) 
z{l - z){z - y){x - z) 
1 

E - qi -k2 - LU-{ki - qi) 
1 



E-Vi 



k2 — CO (fci 

1 



X ■ 



+ {1^2}. 



92 (g 



Pi] 



(51) 



The symbol {1^2} means that all labels 1 are replaced with 2 and vice versa, as well as 
replacing the Bjorken variables x, y, and z with 1 — x, 1 — y, and 1 — z. This is a way 
of explicitly stating the symmetry of the potential under exchange of particles 1 and 2. A 
detailed discussion of the evaluation of the loop integral in Eq. (^ip is given in Appendix 0. 
It is straightforward to calculate the other parts of the TBE potential, 



Vtbe:Sx(^; ki;pi] 



'eJ J 2(27r)3 



Ojx - z)e{z - y) 
z{x — z){l + z — X — y){z — y) 
1 

E -q{ -k2 - u-{ki - qi) 
1 



X 



E -p^ -k2 - id (fci - qi) - iu~{qi - Pi] 
1 



X ■ 



E -Pi -q2-uj {ki- qi] 
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(a) VtbeiSB 


k2 P2 k2 P2 

/ / \ \ 

ki-qi//qi-pi + fc2-g2\\^2-p2 

— ^-^ 

ki qi Pi h qi Pi 


(b) VtBE:SX 


h P2 k2 J?2^ P2 

/h-q, + 
y^qi-pi v\^2-?2 
ki qi Pi h q. Pi 


(c) VtBE:TX 


h , ^2 ^ P2 ^2 , ^2 ^ p2 

y 1^1 Qi _|_ y 

i ^ 

ki qi Pi fci q. Pi 


(d) VtBE:WX 


^2 /2^ P2 k2 , ^2 ^ P2 

\/ \ / 

ki - qi /^\ Pi - qi + ^2 - <?2\ /P2 - ?2 

/ \ Y 
— ^ ^ 

ki qi Pi fci Pi 


(e) VtbeiZX 


k2 P2 k2 ^ ^2 ^ P2 
^^^X^^ _^ ^2 - 92 \ / P2 - ?2 

— — ^i-k— ^ \ 

fci gi Pi ki Pi 



FIG. 2. The components of the TBE potential, (a) the stretched-box, (b) stretched-crossed, (c) 
T-crossed, (d) wide-crossed, and (e) Z-crossed diagrams. Here, qi = ki+pi—qi and q2 = k2+P2—q2- 



+ {1^2}, (52) 

where we have denoted the hght- front energy of particle 2 by §2^, given by 

q^ = e~{P + qi-pi- ki) (53) 

The rest of the TBE potential is given by 
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VTBE:Tx{E;ki;pi) = ( — 



MV 
'EJ 



Vtbe:Zx(-E; ki]pi) 



2(27r)2 



1 ^J{x - z)0{l + z-x- y)e{y - z) 
z{x — z){l + z — X — y){y — z) 



1 



X 



X 



E-Qi -k^ -uj (fci - qi] 
1 



E - qi -q2 {ki-qi)-uj {pi - Qi) 

1 



X- 



r d^qs_ 
'EJ J 2(27r)2 



+ {1^2} 



1 ^ 6'(s - z)6'(l + z - X - y)6{y - z) 
z{x — z){\ + z — X — y){y — z) 



E -q{ -k2 - uj-{ki - qi) 
1 



(55) 



X 



E-qi -q^-uj {ki-qi)-uj {pi - qi] 
1 



X- 



'eJ J 2(27r)2 



+ {1^2} 

6'(x + ?/ — 1 — 2;) 



(56) 



z{x — z){x + y — 1 — z){y — z) 
1 



X 



X 



E -qi -UJ (fci - qi) 

1 



E - gi - -P2 -e-{pi + ki- P -qi- 
1 



X- 



-P2 {Pi -Qi) 
+ {1^2}. 



(57) 



The loop integrals in the expressions for the TBE potentials and the azimuthal-angle aver- 
aging are discussed in Appendix 0. 



C. TBE:SB Potential: Connection to the ladder Bethe-Salpeter equation 

It is well known that the full, untruncated Bethe-Salpeter equation can be reduced to the 
full, untruncated Hamiltonian (Schrodinger-type) equation by integration over the energy 
or light-front energy variables. If a truncated kernel is used for the Bethe-Salpeter equation, 
then the physically equivalent Hamiltonian equation will not include all the graphs that 
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the full theory allows. By physically equivalent, we mean that the spectra of the potential 
V should reproduce the spectrum for the states of the Bethe-Salpeter equation, excluding 

|26[] . For an extensive discussion of this in the equal-time 

Lahiff and Afnan HSOl, and for 



the so-called "abnormal" states 
case see, for instance, Klein 



281, Phillips and Wallace |29 



24] 



examples on the light front, Chang and Ma ||31| and Ligterink and Bakker 

In particular, consider the Bethe-Salpeter equation when the ladder kernel is used. The 
physically equivalent light-front potential will not include any graphs where the meson lines 
cross, so to order g"^, the potential is given by (7^VbBE + 5''^Vtbe:SB- Therefore, by considering 
the TBE:SB truncation, we can test how well the light-front Hamiltonian approach approx- 
imates the full ladder Bethe-Salpeter equation. This is idea discussed more throughly in 



IV. NON-PERTURBATIVE POTENTIALS 

The potentials discussed in this section are derived from the OBE field theory potential, 
but additional approximations are made to simplify the expressions. 

A. Symmetrized-mass approximation 



Krautgartner, Pauli and Wolz ^2[ and Trittmann and Pauli studied positronium 
with a large coupling constant in light-front dynamics. The one-photon-exchange potential 
they obtain has a colinear singularity due to the sum of the instantaneous photon exchange 
graph and a gauge-dependent factor from the spin sum. They argue that the singularity 
is not physical, and therefore must be canceled by higher-order terms in the potential. 
The effect of those terms can be simulated by choosing the bound-state energy so that the 
coefficient of the singular term vanishes. They find that this condition is met when the 
light-front energy P~ in the one-photon-exchange potential is replaced with the operator u, 
expressed here in the two-particle basis. 



P ^uj{ki,k2]Pi,P2) =^{pi +P2 +ki +f^2 



(58) 



This is called the symmetrized mass |3^, the average of the total P~ in the initial and 
final states. It is important to note that this approximation affects not only the singular 
term, but also the energy denominator in the rest of the OBE potential. The modified 
denominators simulate the effects of the non-perturbative higher-order terms that are not 
included explicitly in the OBE potential. Potentials obtained with this approximation are 
similar to those given by the unitary transformation method where the potentials 

depend explicitly on the initial- and final-state energies. 

In our model, there are no singularities associated with the OBE graphs because we deal 
only with scalar fields. However, we may use their approximation to obtain a new light-front 
OBE potential that should incorporate some non-perturbative effects. Recalling that the 
only place where P" occurred in Eq. (^) was in the denominator, the E in the denominator 
of the OBE potential is replaced with uj to get 
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K.(fci;pi) = (-) 



e{x - y)/\x - y\ 



\ {-Pi + P2 + ^1 - ^2 ) - ^ {kl- Pi) 

_^ d{y - x)/\y - x\ 

I {+Pi -P2 - K + A;^) - uj-{pi - kl) 



(59) 



/MY E , , 

— 7 X • 60 

\E) - p+) (-p- + p2 + K - (kl. - 



Writing the light-front variables in the denominator in terms of the equal-time variables, as 
prescribed in Eq. (p2|), we find 



likt - pt) (-pr + P2 + K - k,) = -(P - p^)' + ^^(A:° - p^. (61) 

Thus, Eq. (^0]) can be rewritten as 

/ M\ ^ E 

K.(fcET;PET) = (-^) ^3^3-77 — -■ (62) 

This result can also by obtained more directly by considering the first term in Eq. (^6|). 
Recall that the E'^ that appears in the denominator is written as P^P~ in an arbitrary frame, 
so in the symmetrized-mass approximation, the E'^ term is replaced with Euo. This causes 
A term in the denominator of Eq. (^) to vanish, so the equation reduces to Eq. (|62D . Also, 
note that by writing this new potential, we attempt to incorporate physics from higher-order 
graphs than just the OBE graphs. 

The singularity structure of the symmetrized-mass potential is easily analyzed. When 
scattering states are used, in the center-of-momentum frame the total energy of the state is 
E = 2k^ = 2p^, so the relations in Eq. ( |32D become 

k^ = k^±k\ (63a) 

kf = k^^k\ (63b) 

Using these relations, the symmetrized mass is uj = E. Thus, for scattering states, this 
potential is same as the OBE scattering potential and the singularity structure is the same. 

B. Instantaneous and Retarded approximations 

For our bound states, /c|rp <C M^, so that the p^(fc° — P°)^ will be much smaller 
than the other terms in the denominator of Eq. (0). Therefore, we may approximate the 
symmetrized-mass potential by the instantaneous potential, 

/M\^ -E 

Vlnst{kET;PET) = [-^] 77 T^- ^- (64) 

V ^ / (^ET - Pet)^ + K 
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Alternatively, we may also argue that since the energy difference term is small, we can also 
approximate \{j by the retarded potential, 

VRetikET] Pet) = ^ — IfTT2 7Z Zr^2 (^5) 



(66) 



My E 
EJ (A;et -Pet)^ - /i^' 



where /cet and pet represent four-vectors, defined by the equal-time three-vectors and 
the condition that A;|t = Pet — These potentials resemble the three-dimensional 

Blankenbecler-Sugar or Gross p6| quasi-potentials. 



Both of these approximations are reasonable if the energy difference between the initial 
and final states is small, which is valid for lightly-bound states. The instantaneous poten- 
tial is a better approximation of the symmetrized-mass potential, since if we expand the 
symmetrized-mass potential to second-order in perturbation theory about = we get 
Vuj = ^nst- Also, note that these potentials are explicitly rotationally invariant in terms of 
our equal-time parameterization, which provides significant computational advantages. 



C. Three-dimensional reduction of the Bethe-Salpeter equation 

We now consider a non-perturbative approximation used by Wallace and Mandelzweig 
1 37,38]. The basic idea is to first make an approximation of the Bethe-Salpeter equation, 



then reduce that modified Bethe-Salpeter equation to the physically equivalent Hamiltonian 
equation. This approach was used by Phillips and Wallace |^ for the model we use. 



however, they obtained an equal-time Hamiltonian, while we seek a light-front Hamiltonian. 
Before we do this, we first review the basic mechanics of the three-dimensional reduction. 



as presented in Sales et al. ||2^ and specialized to our particular case. We postpone the 
discussion of the approximation until section [IV D| . 

The Bethe-Salpeter equation can be written in matrix form as F = KGqT, or explicitly 
in function form in the momentum basis as 

nki;P) = J ^^K{h,p,;P)Go{pi;P)T{p,;P). (67) 



i2nY 

In these equations, F is the four- dimensional vertex function, K is the four-dimensional 
kernel, and Go is the two-particle four- dimensional Green's function. The momenta are P, 
the total four-momentum, andpi, the four-momentum of particle 1. Particle 2's momentum 
is implicitly P — pi. The four-dimensional Green's function is given by 

Go{ki; P) = td{ki)d{P - ki), (68) 

where d is the one-particle Green's function. On the light front, d can be written as 

1 




- Sign(p+)e-(p)' 
where the light-front energy e~ is given by 



(69) 
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The real part of is a positive definite quantity, and 1] is positive infinitesimal. 

as 



The Bethe-Salpeter equation can be rewritten |39 



r = WGqT, (71) 
where Go is an auxiliary Green's function, and W is defined by 

W = K + K{Go-Go)W. (72) 
The advantage of this rearrangement is that we are free to choose the form of the auxiliary 



Green's function, On. The choice of Gn advocated in Ref. 23 



IS 



Go{h,p^;P)=Goih;P) ). ^,^ Go{pi;P), (73) 



where 



^foih;P). (74) 

There is an extra factor of 2 in the denominator of Eq. ( [74|) when compared to the equal-time 
formalism. This is due to the Jacobian of the light-front coordinates. 

Using the definition of Gq given in Eq. (0), we can integrate the modified Bethe-Salpeter 
equation, in Eq. (^), over the light-front energy to get 

7(fci;P) = / ^-^^^wik„p,;P)go{pi;PHpuP), (75) 



where 



wik,,p,;P) ^ (GoWG,)ik,,p,;P)— (76) 

9o{ki;P) 9o{Pi;P) 

7(fci; P) ^ —-— / ^Go(fci; P)m; P). (77) 
go{ki;P) J 2(27r) 

The functional (/) is defined by its action on an arbitrary function f{ki,pi), where ki and 
Pi are four-vectors, as 

We proceed by calculating the specific form of go, 

^) = 2 p--fc--fc- ' 
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where k~ = e~{ki). As (70 is a three-dimensional quantity, it is clear that k~ is not the 
independent minus component of a momentum four- vector. With this expression for Qq, we 
can specialize Eq. ( |75D to the center-of-momentum frame and obtain 

(E-k^-k^)^{k,;E)= [d'p^,^ [\pt ^(^y^^'^) J {pi;E), (80) 
^ ^ 7 ^ 2{27ry^ktHptpt 

where 

^(fc,;E) = M^^(fc,;E). (81) 

By comparing this equation to Eq. (|29|), we find that the full light-front two-nucleon effective 
potential that corresponds to the kernel K, after suppressing the coefficient given in Eq. (^Tj), 
is 

Viki,p^;E) = ^wiki,pi;E). (82) 

Thus, we can calculate light-front potentials directly from the Bethe-Salpeter equation using 
this method. 

The potential V can be expanded in powers of the coupling constant, as done in LFTOPT. 
We find that when the auxiliary Green's function given in Eq. (^) and the OBE kernel are 
used, the lowest order parts of the potential (as calculated in [^) are the same as our 
OBE and TBE:SB potentials. Thus, we conclude that this method produces the physically 
equivalent Hamiltonian theory to the Bethe-Salpeter equation being used. We will use this 
in the next section to derive a Hamiltonian potential for a situation where LFTOPT cannot 
be used. 



D. The modified-Green's-function approach 

Now that the technology for the three-dimensional reduction has been reviewed, we 
derive an approximate kernel for the Bethe-Salpeter equation. We will follow the approach 



of Phillips and Wallace and works cited therein. The idea is to start with the Bethe- 
Salpeter equation where the kernel is truncated to only include ladder (one-boson-exchange) 
and crossed (two-boson-exchange) parts, 

r = (-ft'ladder + -f^cross)Gor. (83) 

An uncrossed approximation is used where the crossed part of the kernel is approximated by 
-Across ~ -^ladderCc-f^iaddcr- Our job is to find a Valid modified Green's function, Gc- Using 
this uncrossed approximation, 

r ~ (-ft'ladder + -^ladderGc-f^laddcr)G'or. (84) 



One can attempt to rewrite Eq. ([83| ) as an equation linear in A'ladder? to obtain the 
modified-Green's-function Bethe-Salpeter equation. 
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■TmCF -^ladder 

By iterating this integral equation for Tc, we obtain 



MGF 



oo 



-f^ladder + / K] 



ladder 



ladder J 



n=l 



GoT 



MGF- 



(85) 



(86) 



The part of Eq. (pSD that plays the role of the kernel includes the uncrossed approximation 
of the original kernel -ft'iadder + -^ladder Cc-K^iadder as well as many more terms. We note that 
the higher-order terms approximate some of the higher-order terms that should be included 
in the full kernel, such as three-boson-exchange diagrams where several meson lines cross. 
However, this approach undercounts the higher-order terms which it approximates, and 
also leaves out some terms completely. Therefore, this new Bethe-Salpeter equation will 
give results that are closer to the full solution than Eq. ([83|), but will not give the exact 
solution. The articles by Wallace and Mandelzweig p7||38| demonstrate that this approach, 
by effectively summing an infinite set of interactions, gives the correct one-body limit, which 
is something that the usual Bethe-Salpeter equation with a truncated kernel cannot do. 

The modified Bethe-Salpeter equation in Eq. ( pSf ) is reduced to a Hamiltonian equa- 
tion via the technique discussed in the previous section. The equal-time Hamiltonian has 
been derived by Phillips and Wallace They found that this modified-Green's-function 
approach gave a spectra that lies closer to the full ground-state spectra than the other ap- 
proximations they considered. We will use the light-front reduction to obtain the light-front 
potential for the Hamiltonian equation physically equivalent to Eq. (^). 

To clearly see what role Gc plays, we compare the crossed and uncrossed Feynman 
graphs in Fig. 13. Using the Feynman rules. 



-^crossed ^ 



(27r)4(A;i-gi)2-/i2 
dSi 1 



i^uncros.ed OC j ^^^^^ _ ^^^^ _ 



d{qi)d{q2) 
d{qi)d{q2) 



1 



1 



(87) 

(88) 



where d is the one-particle propagator given in Eq. (pOD, and q2 = —q2 +P2 + ^2- The only 
difference between these two graphs is that the crossed one has d{q2) while the uncrossed 
one has d{q2)- 



Crossed 



?2 



P2 



V 



h - Qi/ - qi 

ki ^ ^» Pi 

qi 



Uncrossed 



92 



P2 



h- qi ! 



Pi - qi 



ki 



I 



qi 



pi 



FIG. 3. The crossed and uncrossed Feynman graphs, where §2 = —Q2 +^2 + ^2- 
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We want an approximate one-particle Green's function dc that only depends on qi and 
P, so that 



dc{qi;P)-d{q2). (89) 

Substitution of dc{qi',P) for d{q2) in the crossed graph causes the graph to become un- 
crossed. The penalty for this simplification is that a modified Green's function propagates 
in the intermediate state, namely id{qi)dc{qi', P)- It is important that this approximation 
is invariant under relabeling particle labels, so we explicitly symmetrize by defining 

Gc='- [rf(gi)rfc(gi; P) + d{q2)dc{q2; P)] (90) 

How should we approximate dc*? Since we are interested in obtaining a bound state, 
a low-energy approximation is chosen. Specializing to the center-of-momentum frame in 
this limit, the external momenta are half the total momentum, so p2 = k2 = P/2 and 
q2 = P — q2 = Qi- This approximation is similar to the one used by Phillips and Wallace. 
Thus, we define dc{qi', P) = d{qi) so 

Gc{k^, P) = G,{k,- P) + G2{k,; P), (91) 

where we define Gi and G2 by 

G,ik^;P)='-dik,y (92) 
G2{k^;P)='-d{P-k^y. (93) 

This approximation for Gc is valid for this model for the energy range we study, as discussed 
in Appendix |^. 

We can write the modified-Green's-function Bethe-Salpeter equation as 

(94) 

where 

Go = Go + Gi + G2. (95) 
This is used as the starting point of three-dimensional reduction discussed in Section [IV C 



where Gq is considered as the Green's function. Before doing the reduction, note that the 
two poles of Gi and G2 lie in the same half plane for each function, so 

/dk^ ~ 
^Go(fci,;P) (96) 

= go{k,,P) (97) 

Proceeding with the three-dimensional reduction of Eq. ( p^ ) in the center-of-momentum 
frame, we obtain 



23 



^ ^ J Jo 2{27rfJktktptpt 



where 



^(kuPuP) = A G,WG,){k,,p,-P)— (99) 

9o{ki;P) 9o{Pi;P) 

~ 1 r dk^ ~ 

tPik,; E) = ^_ / — i-G'o(A;i; P)TuGF{ku P), (100) 



fktkt^ 2(27r 
and the modified kernel W is given by 

W = -ftTiadder + -^ladder 



Go — Gq 



W (101) 



Go{k,,p,;P)=Go{k,-P) lJ_i^Go(pi;P). (102) 

go{ki;P) 

It is a feature of the hght front that Qq = Qo, so that the uncrossed approximation only affects 
the potential, and Eq. (|98|) has the same form as Eq. (ffSl). In the equal-time calculation ||2^ 



9o 7^ fl'O) so the approximation changes both the Green's function as well as the potential. 

We now expand W in powers of the coupling constant, and keep only the lowest order 
term, -ft'iadder- According to Eq. (|52D, the light-front potential that corresponds to this 
truncation of the kernel is the modified-Green's-function (MGF) potential Vmgf, 

VMGF{kl,Pi;P) = ^ ,^ p, (Goi^iaddergo)(fci,pi;P) / ■ (103) 
Ego{ki;P) 9o[Pi]P) 

The one-boson-exchange kernel -ft'iadder is given by the Feynman diagram, so 



Kladdev{kl,Pi; P) 



{ki - piY - ij? + 17] 



\kt -pt) {ki -p^)-Sign{kt -pt)uj-{ki-pi] 
where the light-front energy of the meson is given by 



(104) 



uJ-iq) = ^^-j-^-iV. (105) 

By examining the locations of all the poles in the k~ integrals for Vmgf, we find the 
integrals are no n- vanishing only when both x and y are between and 1. The sign functions 
in the denominator of -R'ladder naturally divide Vmgf into two parts, one for x < ?/ and the 
other for x > y. The integrals in Vmgf are straightforward, but quite lengthy and tedious. 
Therefore, we show only the final answer, 

^^9{x-y)( 1 , Np^2i + Nk,u , Np^2iNk,i2 



W(fci,Pi;P)= (— )" 



\x-y\ \Di 2Df 2Di J 

\y-x\ \D2 2Dl 2DI , 



(106) 
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where 



Nk,u = %iE - - k^) (107) 
K2 

Nk,2i = %{E - k^ - k^) (108) 

D,=E-p^ -k^ -uj-{ki-p2) (109) 
D2 = E-k^ -u;-{pi-k2). (110) 

The expressions for A'p 12 and Np^2i are obtained by replacing k with p in N^^u and A^^^ 21- 

What is the physical interpretation of this modified-Green's-function potential? The 
first term multiplying each 6 function gives the OBE potential we derived before from the 
perturbation theory. There the D in the denominators corresponds to one meson exchange. 
The second and third terms multiplying the 6 functions, with and in the denom- 
inators appear to be effective two- and three-meson-exchange terms. Since time-ordered 
perturbation theory does not apply to the modified Bethe-Salpeter equation that we use, 
the exact nature of these terms is not easy to understand. However, it is clear that these 
terms increase the strength of the potential, and should mimic the higher-order diagrams 
that are not being included explicitly. 

The only dependence on the direction of the perpendicular components of k and p comes 
from the D's. This allows the azimuthal- angle integration of Vmgf to be done easily, as 
shown in Appendix 0. 



V. RESULTS 



For our numerical work, we pick the meson mass to be 0.15 times that of the nucleon, so 
fi = 0.15M. This is chosen so that our ground state can be considered a toy model of deu- 
terium, and also to facilitate comparison with the results of Nieuwenhuis and Tjon and 
Phillips and Afnan Nieuwenhuis and Tjon used the Feynman-Schwinger representation 
(FSR) of the two-particle Green's function in the quenched approximation without the 
mass and vertex renormalization terms . Their result is to be considered the full solution 
that the Bethe-Salpeter and Hamiltonian equations approximate. For the Bethe-Salpeter 
equation, computation of the bound-state energies for models similar to ours have been done 
for the ladder |^l| and ladder plus crossed [jl6| kernels over 30 years ago. More recent results 
are found in 



where the solutions are compared to those given by the FSR approach. 
Now consider how the light-front Hamiltonian approach fits in with the other approaches. 
As discussed in section 111 G , different light-front potentials can be derived from Bethe- 
Salpeter equations with different kernels. We have mentioned that the OBE+TBE:SB poten- 
tial should approximate the ladder Bethe-Salpeter equation, and similarly the OBE+TBE 
potential should approximate the Bethe-Salpeter equation when the ladder plus crossed 
kernel is used. The best that these truncated Hamiltonians can do is approximate their 
respective Bethe-Salpeter equations. 

With this in mind, we evaluate the coupling constant versus bound-state energy curves 
(which we will call the spectrum) for the Hamiltonian equation with the OBE potential, 
the OBE+TBE: SB potential, and the OBE+TBE potential. For the range of values we use 
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here, we find numerical errors in the value of g"^ are less than 2%. Our results (without error 
bars) are plotted along with the results obtained with the ladder BSE [40|, and ladder plus 



crossed BSE |^3| in Fig. ^. We note that the OBE+TBE:SB potential agrees well with the 
ladder BSE, and the OBE+TBE potential agrees with the ladder plus crossed BSE. This is 
the best that a Hamiltonian can do, so this result is interpreted as evidence that, in general, 
the higher-order diagrams are very small for the ground state on the light front. 
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FIG. 4. We display here three light-front calculations using the OBE, OBE-|-TBE:SB, and 
OBE-I-TBE potentials. The spectra for the Bethe-Salpeter equation for the ladder and ladder plus 
crossed kernels are also plotted. The curves for the ladder (ladder plus crossed) Bethe-Salpeter 
equation and the light-front OBE-|-TBE:SB (OBE-I-TBE) potentials are very close to each other, 
almost indistinguishable in this figure. E is the energy of the ground state of two nucleons, and 
M is the mass of the nucleons. The meson mass is ^ = 0.15M. In terms of the binding energy i?, 
E = 2M -B. 



If all one wanted was a way to approximate the spectra for Bethe-Salpeter equations with 
different kernels, one could just use the truncated potentials that the LFTOPT provide. 
However, the true goal is to approximate the spectra for the full ground state, which in this 
model is given by the FSR approach We expect that the non-perturbative potentials 
should give a better approximation of the full solution than the perturbative potentials, since 
the non-perturbative potentials attempt to incorporate physics from higher-order diagrams, 
although this is not immediately clear by looking at the forms of the potentials used. We plot 
the results for all of the light-front potentials described in this paper, the three truncated 
potentials (OBE, OBE-|-TBE:SB, and OBE-f TBE) and the four non-perturbative potentials 
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(symmetrized mass, retarded, instantaneous, and modified Green's function) along witli tlie 
results for the full theory in Fig. ^. For deeply-bound states, there is considerable dis- 
agreement between the perturbative results and the full results, while the non-perturbative 
results do better, with the modified-Green's-function (MGF) potential achieving the closest 
agreement. For lightly-bound states, the results for all of the potentials appear converge to 
each other, close to the full result. 
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FIG. 5. The spectra for the seven light-front potentials (QBE, OBE+TBE:SB, OBE-hTBE, 
retarded, instantaneous, symmetrized mass, modified-Green's-function) are plotted here along with 
the spectra for the full solution (FSR). The instantaneous curve lies very close to the symmetrized 
mass curve, hence the same line style is used for both. 



For the modified-Green's-function potential, only the first term of the expansion in was 
kept. In principle, higher-order terms could be calculated. However, since there was fairly 
good agreement between the OBE potential and the ladder Bethe-Salpeter equation, it is 
expected that the MGF potential will give results that are close to the ladder Bethe-Salpeter 
equation using the modified-Green's-function. 



VI. CONCLUSIONS 

In this paper, we consider the ground state of a massive Wick-Cutkosky model using a 
Hamiltonian derived from light-front dynamics. We examine three different truncations of 
the effective potential derived from the perturbative field theory, and four approximations 
that attempt to incorporate non-perturbative physics. For each of these potentials, we 
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calculate the coupling constant that gives the ground-state for a given bound-state energy 
(the spectra), and compare to the spectra from different approaches found in the literature. 
We find fairly good agreement between all the methods for lightly-bound systems. 

For the full range of binding energies studied, the results for calculations including one- 
and two-boson-exchange potentials agree with the Bethe-Salpeter equation using the physi- 
cally equivalent truncation of the kernel within the numerical errors. (This is a consequence 
of examining the ground state. For the excited states more higher-order light-front time- 
ordered graphs are required to get the same level of agreement 0.) The agreement for the 
case with the stretched box diagrams (OBE-|-TBE:SB) has been shown previously by Sales 
et al. the result with the crossed-box contribution (OBE+TBE) is new. This excellent 
agreement with the Bethe-Salpeter results has an undesirable consequence: The BSE results 
are known to be a poor approximation of the full solution for deeply bound systems, so the 
truncated Hamiltonian approach cannot provide a good approximation to the full solution 
in that regime. 

The non-perturbative potentials based on physical considerations give a better approxi- 
mation of the full solution than the potentials obtained from LFTOPT. For all binding ener- 
gies, the modified-Green's-function potential achieves the closest agreement with the full so- 
lution of all the potentials considered here. However, there is still considerable disagreement 
between approximate potentials and the full result for deeply bound states. We interpret 
this as an indication that the approximations, while incorporating some non-perturbative 
physics, do not go far enough. In the weakly-bound regime, which is of relevance for deuteron 
calculations, the spectra for all of the potentials are close together, indicating that light-front 
dynamics provides a good description of lightly-bound systems. 
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APPENDIX A: NOTATION, CONVENTIONS, AND USEFUL RELATIONS 



ACKNOWLEDGMENTS 



This is patterned after the review by Harindranath [Q. For a general four- vector a, we 
define the light-front variables 




so the 4-vector can be denoted 



a = [a ,a ,a±). 



(A3) 



Using this, we find that the scalar product is 




(A4) 



28 



This defines g^^, with = = 1/2, gn = g22 = —1, and all other elements of g 
vanish. The elements of g'^'^ are obtained from the condition that g'^'^ is the inverse of (yf^j^, 
so g'^^gpx = 5". Its elements are the same as those of g^y, except for g~^ = g^~ = 2. Thus, 

= 2a^. (A5) 

and the partial derivatives are similarly given by 

a± = 2d^ = 2-^. (A6) 

Moving to the physical consequences of this coordinate system, the commutation rela- 
tions \p^, x"] = ig^^ yields 

[p^, x^] = 2i (A7) 
= (A8) 

with the other commutators equal to zero. Thus, x\ is canonically conjugate to p^, and x^ 
is conjugate to . In light-front dynamics, x"*" plays the role of time (the light-front time), 
so p~ is the light-front energy and the light-front Hamiltonian is given by P~ . 

Particles have the light-front energy defined by the on-shell constraint k"^ = m?. This 
implies that the light-front energy is 

k- = !^f±^i. (A9) 

The free components of the momentum can be written as the light-front three-vector A;, 
denoted by 

k = {k+,ki_). (AlO) 



APPENDIX B: CONVERSION TO MATRIX FORM 

To solve for the bound-state wavefunction numerically, the light-front Schrodinger equa- 
tion given in Eq. (37) must be discretized and cast in matrix form. The equation is first 
symmetrized to get 

~L "^^^^i ^W(^ET,^fc;PET,^p)^f (PET,^p) =V^f (fcET,^fc), (Bl) 

where 



V^ikET, Ok] Pet, Op) = B{kET)A{kET, Ok)V+{kET, Ok^PET, Op)A{pET, Op)B{pet) 

V^f (fcET,4) = B{kET)-'A{kET,Ok)ij''^{kET,Ok) 



A{kET, Ok) 

B{kET) 



\ 



sin Oh 

W) 



0\2 ■ 



(B2) 
(B3) 

(B4) 
(B5) 
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Before discretizing the integrals, note that 



dpfip) = a['du (fiau) + ^^^] . (B6) 
JO V 



Using this trick, the pet integral in Eq. can be written as an integral over a finite range. 
Since we are concerned with a bound state, the wavefunction is exponentially damped for 
large momenta, and the second term of Eq. ( P6|) converges as u approaches zero. 

All the integrals in Eq. ^T] then are over a finite range, and can be discretized using 
Gauss-Legendre quadrature. The specific routines for the quadrature are given by Numerical 
Recipes in C [Q]. This conversion gives a matrix equation that approximates the original 
Eq. 

-Vs'-i9iE),E)^l;f = ^pf, (B7) 

where the explicit dependence of Vg' on the symmetrized potential on the binding energy 
E and the coupling constant g is shown. This equation must be solved self-consistently for 
the spectrum g{E). 

The approach we use is to first solve for the spectrum for the OBE potential. The 
eigenvalue equation 

^5toBE(^)^f = , (B8) 

where 

a = ~] (B9) 
The ground-state wavefunction is the eigenvector that corresponds to the smallest eigenvalue 



a. We calculate the wavefunction and the smallest coupling constant using EISPACK |^ 
routines for a range of energies to map out the spectrum. 

Using the coupling constant for the OBE potential as a starting point, we can use Eq. |B^ 
for higher-order potentials that include N meson exchanges. For a given energy, the coupling 
constant g{E) is initially chosen as goBE{E), then we solve 



ijf = P^f, (BIO) 



as an eigenvalue equation for p. The coupling constant g{E) is varied until the the lowest 
eigenvalue is /3 = — 1, at which point g{E) is the correct value of the spectrum corresponding 
to the ground-state wavefunction ^ipg^. 



APPENDIX C: AZIMUTHAL- ANGLE INTEGRATION OF THE OBE AND MGF 

POTENTIALS 

In this section, we evaluate the azimuthal-angle integration of the OBE potential in 
Eq. ( ^5|) and the first term in MGF potential in Eq. ( |106| ), using the prescription for 
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azimuthal-angle integration given in Eq. (0). One of the integrals is easily done since 
since the potential is independent of the azimuthal angle between the two perpendicular 
momenta, so 



V{k~^,k±;p'^,p±) oc 
where 



e{x - y) 



2n 



Ai + B cos ( 



+ 0{y - x) 



2n 



A2 + B cos ( 



Ai = {kt - vt) {E-Vi- -fi^-pl-kl 
A, = {pt - k+KE - k^ - p,) - /x^ - pI - kl 
B = 2k±p±. 

The integrals in Eq. are easily done to give, since the A's are negative, 

r'^^ dd) -27r 



/o A + B cos (j) ^A^ - fi2 
Using this, the azimuthal-angle-averaged OBE potential is given by 



VoBEik^,kr,P^,P±) = -2n(—j E 



[x-y) e{y - x) 



Al-B^ 



Al-B^ 



(CI) 



(C2) 
(C3) 
(C4) 



(C5) 



(C6) 



It is straightforward to rewrite this equation for the potential in terms of the equal-time 
coordinates. 

When the other terms in the MGF potential are azimuthal-angle averaged, integrations 
similar to the one given in Eq. are encountered, with the denominator squared or cubed. 
We note that 



27r 



-27r 



A 



(A + 5 cos 0)2 _ 52 ^2 _ 

2^ dS -27r 2^2 + 52 



/O (A + 5 COS 0)3 _52 (A2 -52)2^ 

SO the azimuthal-angle-averaged MGF potential is given by 



(C7) 
(C8) 



l^MGF(fci,pi;P) = -27r(^— 



x-y) ^ A/p,2i + Nk,i2 A/p,2iA^fc,i2 



Al-B^\ 

djy - x) 
Al - fi2 



2Di 



,3 



ID 



2,2 



ID 



2,3 



(C9) 



where 



A,2 = 
A,3 = 



A2-g2 
A 

(^2 - ^2)2 
2^2 + 52 ^ 



(CIO) 

(Cll) 



and i = 1,2. 
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APPENDIX D: AZIMUTHAL- ANGLE INTEGRATION AND LOOP 
INTEGRATION OF THE TBE POTENTIALS 



As in the previous section, we want the azimuthal-angle integrals of the TBE potentials 
given in Eqs. (|5Tl-p7|). For these potentials, there is also a loop integral that has to be done. 
We start by analyzing the equations schematically. Each of the terms in the TBE potentials 
can be written in the following form. 



^TBE(A;+,fcx;p+,P±) = |^'rfzJ(A;+,g+,p+)/(A;+,A:^,g+,g^,p+,p^) (Dl) 

/(/c+,/c_L,g"^,g±,p"^,p±) 



^ Ai+ Bi COS(f)q 

1 

X 



A2 + B2 COS (j)q + C2 COS 0p + D2 COs(0p - 0q) 

X \ ^ (D2) 

A3 + S3 COS 0g + C3 COS 0p + 1^3 COs(0p - 0g) ' 

where the A's, S's, C's, and i^'s may have dependence on A;"*", k^, p~^, p±, = zE, and q±; 
they are independent of the azimuthal angles. These functions can be easily determined for 
each potential by examining the forms of the original equations. The rotational invariance 
of the potential about the three-axis allows the (pk integration to be done trivially. 

In the integrand of /, only the last two terms depend on (pp. To emphasize this, we write 

I{k^,k±,q^,q±,p^,p±) = d(pq — 03 

Jo {Ai + Bi cos (pq) {A2 + B2 cos (pq) [A-^ + B^ cos (pq) 

l2ik^,k^,q,p\p^) = r -— , ^, . J^^' , ^, . , , (D4) 

JO (1 + 02 COS (pr, + 02 sm 0„ ) ( 1 + 03 COS (pp + O3 sm (pp) 



where, for i = 2, 3, 



Ci + Di cos (pq 
Ai + Bi cos (pq 
Di sin (pq 



Ai + Bi cos (pq 
The integral in I2 is evaluated to obtain 

27r 



(D5) 
(D6) 



l2ik~^,k±,q,p'^,p±) 



(02 - 03)2 + (62 - 63)2 - (a2&3 - 03^2)^ 



^ I 02(^2 - as) + ^2(^2 - ^3) _|_ ^3(^3 - ^2) + ^3(^3 - h] j 



The remaining three-dimensional loop integral in Vtbe on q is done using numeric tech- 
niques. The trick introduced in Appendix ^ to convert the semi-infinite q± integration into 
an integration on a compact range. Before doing the z integral, the range of integration is 
limited by using the step functions. Gauss- Legendre quadrature, given by Numerical Recipes 
in C , is used to evaluate all the integrals. 
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Since each of the parts of the full TBE potential (TBE:SB, TBE:SX, . . . ) should be 
hermitian and invariant under interchange of particle 1 and 2, these invariances can be 
used as a self-consistency check. Each matrix element is calculated twice, first by using the 
straightforward approach, then particle labels 1 and 2 are interchanged and it is calculated 
again. The results are compared, and if they differ by an unacceptable amount, the number of 
quadrature points is increased and the element is recalculated. In order to get the numerical 
accuracy of the potentials correct to within 1%, we start with ten points for the q± integral, 
six points for the (j)q integral, and three points for the z integral, resulting in a three- 
dimensional integral using 180 points. 



APPENDIX E: CHECK OF THE UNCROSSED APPROXIMATION 



In this section, we want to check that how well the approximation 



KlsAdevGcKlaAder , 



(El) 



works. Since we are using a Hamiltonian theory and are interested in the potentials, we 
compare the potentials defined by 



VtbeiX — -^do ^ {GoKcrossGo)go ^ 

VtBE:UX = -^90 ^ (GoKis^ddcrGcKiadderGo) Qq ^. 



(E2) 
(E3) 



The notation used here is defined in sections 11 V CI and 11 VD 



The TBE crossed potential (TBE:X) can be written as 

VtBEiX = ^BE:SX + ^BE:TX + VtBEiWX + ^BE:ZX5 



(E4) 



where the potentials on the right-hand side are defined in section pi B| . Calculation of the 
TBE approximate uncrossed potential (TBE:UX) is straightforward, but tedious. We find 
that 



Vtbe:Ux(-E; ki,pi] 



V'tbe:Uxi {E;ki,pi) + KtbeiUxi Pi , fci ) 

+ ^TBE:UX2(^; + ^TBE:UX2 Pi , ^l) 



(E5) 



where 



VtBE:UXi(^; = (^ — ^ 



\EJ J 2(27r 



^ ^J{x - z)e{z -y) 
{x — z)z'^{z — y) 



1 



X ■ 



E -pi - k2 -uj {qi -pi] 
( 1 



uj {ki - qi) 

\ 2' 



\E-qi 
+ {1^2} 



UJ (fci - qi) 



(E6) 
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EJ J 2(27r) 



^^ 9{x - z)e{y - z) 
10 {x — z)z'^{y — z) 
1 

X — 

E-qi -P2 - ^ [Pi- qi 



\E - Qi -UJ (fci - qi) 

+ {1 ^ 2}. 



(E7) 



Now consider the azimuthal-angle and loop integrals for these potentials. The approach 
used is similar to that of section 0. Analyzing the potentials reveals that each of the terms 
can be written in the following schematic form, 



VTBE:VX,i{k^,k_i]p+,p^) = / dz J{k+,q+,p'^)Ii{k^,k_i,q^,q_L,p+,pjL) 

Jo 2[2tt)'^ Jo 

Ii{k^,k±,q^,q±,p^,p±) = 

Jo 



Ai i + Bi i cos (j)q + Ci^i COS (j) \ A2 + B2 COS ( 



(E8) 
(E9) 



where 



= + 0p, and the A's, -B's, and C's may have dependence on fc"*", A;j_, 
g+ = zE, and g^; they are independent of the azimuthal angles. These factors can be 
easily determined for each potential by examining the forms of the original equations. The 
rotational invariance of the potential about the three-axis allows the 0^ integration to be 
done trivially. The integral is easily done to obtain 



A;_L,g"^,g±,P^,P±) = -27r 



27r 



1 



where 

tti = Ai^i + Bi^i COS05. 
Further simplification is possible for Vtbe;UX2, since for that potential Bi^2 



I{k+, k±,q'^,q±,p'^,p±) 



(27r)2A2 



A12 ~ C^l,2\/ ^2 ~ -^1(^2 



B^) 



(ElO) 



(Ell) 



(E12) 



The techniques discussed in the previous section are used to do the remaining loop integrals. 

The spectra for the OBE+TBE:SB+TBE:UX potential can be calculated and compared 
to the OBE+TBE:SB+TBE:X potential (which is the same as the OBE+TBE potential), the 
modified-Green's-function potential, and the ladder plus crossed Bethe-Salpeter equation. 
The spectra are plotted in Fig. ^ The spectra for the TBE:UX, TBE:X and BSE all lie 
close to each other, which shows that that the uncrossed approximation is valid. 

This shows that the important approximation in the modified-Green's-function approach 
is not the uncrossed approximation, but the addition of the extra interaction terms in 
Eq. ( ^5|) which serve to mimic the higher order interactions. 
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FIG. 6. The spectra are plotted for the potential derived from the TBE truncation of the un- 
crossed approximation (OBE+TBE:SB+TBE:UX, denoted in the figure by TBErlJX), the TBE 
potential (OBE+TBE, denoted by TBE:X), the ladder plus crossed Bethe-Salpeter equation 
(BSE(L-I-X)), and the modified-Green's-function potential (MGF). Except for the MGF spectrum, 
all the curves lie close to each other. 
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